Universal Critical Behavior of the Synchronization Transition in Delayed Chaotic Systems 
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We numerically investigate the critical behavior of the synchronization transition of two unidirectionally 
coupled delayed chaotic systems. We map the problem to a spatially extended system to show that the synchro- 
nization transition in delayed systems exhibits universal critical properties. We find that the synchronization 
transition is absorbing and generically belongs to the universality class of the bounded Kardar-Parisi-Zhang 
equation, as occurs in the case of extended systems. We also argue that directed percolation critical behavior 
may emerge for systems with strong nonlinearities 
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Synchronization of chaotic systems has attracted much in- 
terest in recent years and examples include chemical reac- 
tions, neuronal networks, Josephson junctions, electronic cir- 
cuits, and semiconductor lasers, among others (see Ref. 1 1 ] 
and references therein). More recently, and from a practical 
point of view, this burst of activity in the field is partially due 
to the potential applications in control and secure communica- 
tions. It is expected that an increased complexity of the attrac- 
tor would make it much more difficult to extract the dynamical 
information. In particular, delayed dynamical systems have 
been suggested as the ideal candidates for secure communi- 
cation for several reasons. On the one hand, they are hyper- 
chaotic systems with an arbitrarily large number of positive 
Lyapunov exponents, whose number increases linearly with 
the delay time 0] ■ On the other hand, they may be exper- 
imentally realized in the form of fast communication optical 
systems byusing different types of delayed feedback setups 

Synchronization of two separate delayed chaotic systems 
is achieved by allowing some communication between them. 
The possible schemes are diverse, including variable substitu- 
tion, symmetric feedback, etc. Generally, there exists a critical 
coupling constant k c that separates two different phases. For 
low coupling values, k < k c , there is a disordered phase in 
which each system evolves independently and the time aver- 
age difference between both systems remains finite. In con- 
trast, a synchronized phase appears for k > K c , in which the 
average error tends to zero and memory of the initial differ- 
ence is asymptotically lost. 

Very recent studies have been devoted to investigate im- 
portant aspects of synchronization in delayed dynamical sys- 
tems as, for instance, analytical approximations to estimate 
the synchronization threshold 1 1 1], the robustness of the tran- 
sition to parameter mismatch 1 12], chaos control in lasers with 
feedback 1 13], information flow between drive and response 
systems 1 1 411 - and the effect of a time-dependent delay 11511 . 
However, the mechanism behind the synchronization transi- 
tion in delayed dynamical systems and its relationship with 
those exhibited by other chaotic systems with many degrees 
of freedom is still unknown. 

Remarkably, synchronization also takes place in coupled 



spatially extended systems with many degrees of freedom and 
space-time chaos. In this case, the synchronization transition 
has been shown to be an absorbing nonequilibrium phase tran- 
sition and, accordingly, its critical properties have attracted 
much interest in the past few years fid fnl 111 fll |2(i [ EUl . 
Despite being scalar (i.e., described by only one dynami- 
cal variable), delayed dynamical systems are formally infi- 
nite dimensional dynamical systems and show many aspects 
of space-time chaos, including the formation and propaga- 
tion of structures, defects and spatiotemporal intermittency 
on Him. An interesting question that naturally arises 
is whether the synchronization transition in (scalar) hyper- 
chaotic systems with delayed feedback could also be under- 
stood as a nonequilibrium phase transition, as occurs in (vec- 
torial) extended dynamical systems with space-time chaos. 

In this Letter, we characterize the synchronization transi- 
tion in unidirectionally coupled delayed dynamical systems 
as a nonequilibrium critical phase transition and relate it to 
existing universality classes. We exploit the interpretation 
of delayed dynamical systems as spatially extended systems 
I22I E3I I24I E5I1 to show that the synchronization transition 
in delayed systems exhibits universal properties, which are 
independent of microscopic details of the individual systems 
being coupled. We find that the synchronization transition 
generically belongs to the universality class of the bounded 
Kardar-Parisi-Zhang (BKPZ) equ ation, as occurs in the case 
of extended systems CB 03 El E El . We also argue 
that directed percolation (DP) critical behavior may emerge 
for systems with strong nonlinearities. Our results show that 
the critical properties of the synchronization transition in de- 
layed chaotic systems are identical to those in spatially ex- 
tended systems, despite being the former a scalar system with 
no real spatial structure. 

We consider two identical time-delay systems described by 
two coupled differential-delay equations, the drive (transmit- 
ter) system 



u = —a u + F(u T ), 
and the response (receiver) system 

v = —av + F{v T ) + k{u — v), 



(1) 



(2) 



2 



(a) 

0.07 
0.06 
0.05 
0.04 
0.03 
0.02 
0.01 
0.00 



(b) 



BO .4. 

o 





k= 0.705 




— K= 0.7075 




— K= 0.7025 





0.60 0.64 0.68 0.72 
K 



12 3 4 



FIG. 1 : (Color online) In (a) the average synchronization error in the 
stationary state, w s , is plotted near the synchronization threshold for 
t = 200. Each point corresponds to an average over 600 realizations, 
(b) shows the spatial average of the synchronization error, w(6), for 
r = 2000 and 3 different values of k, each curve being an average 
over 300 realizations. k c = 0.705 is obtained and the corresponding 
slope yields 5 = 1.14. 



where u T = u(t — r) and v T = v(t — r) are the delayed 
variables, r is the time delay, and k is the coupling strength. 

Delayed systems like Eq.Q are used in a variety of appli- 
cations ranging from biology to optics. We have studied in de- 
tail three prototypical models: the Mackey-Glass model 1 26], 
F(p) = b p/(l + p 10 ), (initially introduced to describe regu- 
lation of blood cell production in patients with leukemia), the 
Ikeda differential-delay equation |4], F(p) = b sin(p — po) 
(which appears in the context of optical feedback on a laser- 
beam 1 4] and experimental setups of optical generators of 
chaos in wavelength |8l |9[) and a model with the piecewise 
linear delay-expression given by F(p) — 2p if p < 1/2 and 
F(p) = 2 — 2p if p > 1/2. We have studied the synchroniza- 
tion critical properties by means of computer simulations of 
these three systems and found similar results. 

In order to make apparent the existence of a nonequilibrium 
phase transition we transform the pair of coupled delayed sys- 
tems Eqs.Q and (0 into two coupled spatially extended sys- 
tems. This can be readily done by introducing the coordinate 
transformation, t = x + 9t, where x 6 [0, r] is the space 
variable, while 9 e N is a discrete time variable J22|. Note 
that the time delay r becomes the system size, in such a way 
that the time dependence with the delayed variable is trans- 
formed into an interaction within the horizontal space coor- 
dinate x in the space-time representation. This is a powerful 
representation in which delayed systems can be treated as ex- 
tended systems to identify many features of space-time chaos 

Complete synchronization of the two coupled delayed sys- 
tems, Eqs.Q and (|2j, is achieved if the synchronization error 
u{t) — v(t) vanishes for all times t as t — > oo. In the spatial 
picture we replace the dynamical variables u(t) and v(t) by 
u(x, 9) and v(x, 9), so that the synchronization error is given 
by w(x, 9) — u(x, 9) — v(x, 9), and synchronization occurs 
when w(x, 9) vanishes at all x for 9 — > oo. This is equiva- 
lent to a vanishing spatial average (\w(x, 9)\) x . Note that the 
spatial average in the space-time representation corresponds 
to the average within the delay time r. In contrast, one has 
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FIG. 2: (Color online) Finite-size data in (a) unsealed and (b) 
scaled coordinates for r £ {10, 20, 50, 100, 200, 150, 500, 2000} 
and Kc = 0.705. The collapse is obtained for S = 1.15 and 
z — 1.45. Every curve corresponds to an average over 60 realiza- 
tions. 



(\w(x, 9) \)) x > in the unsynchronized state. This makes the 
average error w(9) = (\w(x, 9)\) x a natural order parameter 
for the transition. Critical properties of the synchronization 
transition can now be studied by analyzing the dependence of 
the order parameter on the coupling strength k. In addition, 
the analysis of the critical behavior for finite time delays can 
be naturally carried out by standard finite-size scaling tech- 
niques. The remaining part of this Letter is devoted to the 
study of these issues. 

Our findings are based upon extensive numerical simu- 
lations of the three time-delay systems introduced above. 
In all our numerical simulations we have used the Adams- 
Bashforth-Moulton predictor-corrector scheme 12711 to inte- 
grate the coupled differential-delay equations and Q. For 
the sake of brevity we focus the discussion of our numerical 
results on the Mackey-Glass model, but we found similar re- 
sults for the Ikeda equation and the piecewise linear system. 
The parameters a = 1 and 6 = 2 are used in all the results we 
are presenting here, and simulation with a time delay varying 
from ten to a few thousand time units have been carried out 
using an integration step of At = 0.01. The region of interest 
here corresponds to delays r ^ 1.7 for which the Mackey- 
Glass model is hyper-chaotic 0]. 

In Fig. la we present our results for the order parameter, 
i.e., the average synchronization error, in the stationary state 
w s (k) = \img^ OQ (\w(x, 9)\) x q for a system of size (de- 
lay) t — 200 as the coupling strength is varied. Inspection 
of Fig. la indicates that the transition is continuous and oc- 
curs roughly around k — 0.7, which is in agreement with 
earlier estimations for the same model parameters II ill . Dy- 
namic critical behavior is studied by calculating the indices 
S and (3 that describe the critical behavior of the order pa- 
rameter near the threshold for synchronization: w{9) ~ 9~ s 
for k = k c and w s (k) ~ \k — k c \@ as the transition is ap- 
proached for k < k c . Studying critical behavior demands to 
obtain a good estimation for the critical threshold which im- 
plies the use of large delays (system sizes). In Fig. lb we plot 
the sub-critical and super-critical behavior of the average syn- 
chronization error for a large system size r = 2000 as the 
transition is approached from below (k = 0.7025) and above 
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FIG. 3: (Color online) Off-critical data in (a) unsealed and (b) scaled 
coordinates for r = 2000 and various values of k £ [0.685, 0.725]. 
The collapse is obtained for k c = 0.705, 8 — 1.05 and V\\ — 1.4. 
Every curve corresponds to an average over 60 realizations. 

(k = 0.7075), respectively. Within our numerical resolution 
we find that the best power-law behavior w(9) ~ 0~ s is ob- 
tained at k c = 0.705 ± 0.002, which gives an estimation for 
the critical exponent 6 — 1.14 ± 0.03. 

Once the critical threshold has been determined, we use 
finite-size scaling at the critical point k c = 0.705 and fit nu- 
merical data to the scaling form 

w{6) = 0- 5 f(6/T z ), (3) 

where the scaling function f(y) ~ const for y <C 1 and 
f(y) ~ y S f° r U ^ 1- This gives us an independent deter- 
mination of 8 and the dynamic exponent z. In Fig. 2a we 
show numerical results for different system sizes at the criti- 
cal point k c and these data are best collapsed in Fig. 2b with 
z = 1.45 ± 0.05 and 8 = 1.15 ± 0.05 (the latter in good 
agreement with our previous estimate in Fig. 1). 

Next we report on off-critical numerical calculations of the 
synchronization error. This allows us to estimate the correla- 
tion length exponent. In Fig. 3a we plot the order parameter 
as the coupling strength is varied close to the synchronization 
threshold for a large system size r = 2000. Near and below 
the transition the characteristic size of synchronized regions 
within the disordered phase is given by the horizontal correla- 
tion length £ and is expected to diverge as £ ~ e~ v± when the 
distance to the critical point tends to zero, e = \n — k c \ —* 0. 
Correspondingly, the characteristic time i3 measuring the du- 
ration of a fluctuation of size £ diverges as d ~ ^ z ~ e _ly " , 
where v\\ = v^z. Off-critical data are then expected to satisfy 
the scaling form w(9, e) = 0~ 5 g(0/i!)), so that numerical data 
in Fig. 3a must collapse according to 

w(0,e) =6- s g{6e v w) (4) 

for the appropriate election of the critical exponents 8 and 
Vn. In Fig. 3b we show a data collapse for the exponents 
5 = 1.05 ± 0.05 and v\\ = 1.4 ± 0.1 where the two branches 
correspond to numerical data for coupling strengths above and 
below critical. Also the index (3 can be obtained from the scal- 
ing behavior Eq. (|4}, for 6 3> d we have w(6 — ► oc, e) ~ e* 3 , 
where /3 = m5 = 1.47. 
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FIG. 4: (Color online) The transverse Lyapunov exponent, Ax, is 
plotted for r = 2000 and various values of « around the synchro- 
nization threshold. Note that Ax changes sign near k c . 



The critical exponents of the synchronization transition in 
delayed systems are to be compared with those observed in 
extended systems. This will allow us to identify the mecha- 
nisms behind the transition in both types of high-dimensional 
dynamical systems. In the context of extended systems, the 
exponential growth rate of the error \w{x, t)\ is known as the 
transverse Lyapunov exponent Aj^ and measures the stability 
of the synchronized solution w(x, t) — 0. Accordingly, stable 
synchronization implies that the transverse Lyapunov expo- 
nent must be negative. For dynamical systems with smooth 
local nonlinearities (like lattices of logistic or tent coupled 
maps), this proves to be a sufficient condition as well. In this 
case, the synchronization transition is found to be generically 
in the universality class of the KPZ equation with a (bound- 
ing) growth-limiting term, the so-called BKPZ universality 
1171 Il8l Il9l l20l bill . Numerical estimates of the critical ex- 
ponents gave Sbkpz = 1.17 ± 0.05, Pbkpz = 1-50 ± 0.05 
and zbkpz — 1-53 ± 0.05 for different models studied in 
the recent literature 1181 Il9t I2C1 \2 ill . On the contrary, in the 
presence of strong and localized nonlinearities (like for in- 
stance for Bernoulli coupled maps), the synchronized phase 
turns out to be unstable even for negative values of Ax O- 
In this case, the transition occurs only when the propagation 
velocity of finite-amplitude perturbations vanishes. The criti- 
cal properties of the transition are then associated with the DP 
universality class flilElll . The fraction of non-synchronized 
sites corresponds to the fraction of active sites in DP. Corre- 
spondingly, the critical exponents are given by 8dp = 0.159, 
(3 DP = 0.277, z DP = 1.581 HE]]]. The DP correlation 
length and time exponents are known to be v± = 1.10 and 
V\\ — 1.73, respectively |28]. 

Our numerical results, 8 — 1.15, (3 — 1.47, z — 1.45 
and m — 1.14, strongly suggest that the synchronization tran- 
sition in delayed chaotic systems generically belongs to the 
BKPZ universality class, as occurs in extended chaotic sys- 
tems. As an additional check, we have measured the trans- 
verse Lyapunov exponent Ax for the coupled delayed sys- 
tems, Eqs. Q and (|2ji, with different coupling strengths and, 
as showed in Fig. 4, we found that the transition takes place 
when Ax becomes negative, as expected for the BKPZ univer- 
sality class. Nevertheless, the nature of the transition can be 
changed to DP behavior in the presence of strong local non- 
linearities akin to what occurs in extended systems. In fact, 
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FIG. 5: The synchronization error for a coupling slightly above the 
synchronization threshold, k > k c , is plotted for r = 500 in the 
cases of (a) smooth nonlinearities and (b) strong local nonlinearities. 
The horizontal width has been chosen slightly larger than r to elimi- 
nate the systematic drift. Compare with Fig.2 of Ref. flHl 

by choosing the nonlinear function F(p) = 2p mod 1 the 
exponent S drops to S = 0.16 ± 0.03 in good agreement with 
DP. The different nature of the transitions in the two cases can 
be seen in Fig. 5, where we show the spatiotemporal evolu- 
tion of the synchronization error \w(x, ff)\ for the (a) smooth 
Mackey-Glass and (b) strongly nonlinear model, respectively, 
just slightly above the transition. 

In conclusion, we have studied the critical properties of the 
synchronization transition in unidirectionally coupled delayed 
chaotic systems. We used a standard coordinate transforma- 
tion to map the (scalar) time-delay system to a spatially ex- 
tended system. This mapping allowed us to study the syn- 
chronization transition as a absorbing nonequilibrium phase 
transition. Comparison of the critical exponents as well as the 
behavior of the transverse Lyapunov exponent lead us to con- 
clude that the synchronization transition in delayed systems 
generically belongs to the BKPZ universality class indepen- 
dently of the specific form of the delay-expression, as long as 
it is a continuous function, just as occurs for synchronization 
of space-time chaos. Finally, our numerical results also indi- 
cate that the existence of discontinuities in the delay nonlinear 
function may change this critical behavior from BKPZ to DP, 
which suggests that the same mechanisms that produce DP 
behavior in coupled extended systems can be invoked in the 
case of delayed chaotic systems. 
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